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ABSTRACT 

Using the Fisher matrix formalism, we quantitatively investigate the constraints on a 10 dimensional 
space of cosmological parameters which may be obtained with future cluster surveys. We explore the 
dependence of the r2 m constraint on both angular coverage and depth of field. We show that in each 
case there is a natural cutoff beyond which the constraints on f2 m do not significantly improve. We 
also investigate the sensitivity of the constraints to changes in our knowledge of the Mass- Temperature 
(M-T) relation by including its normalization and scatter as two of the parameters in the Fisher matrix. 
To make our analysis more realistic, we have added, as priors, the Fisher matrices from hypothetical 
supernova and CMB experiments. We find that X-ray cluster surveys actually help to constrain the M-T 
relation, and explore the implications of this result. 

Subject headings: Galaxies-clusters, cosmology-theory 
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1. INTRODUCTION 

One of the most puzzling mysteries in cosmology today 
is the nature of the dark energy, believed to be driving the 
accelerated expansion of the universe. This constituent is 
predicted to be smooth, except possibly on horizon scales, 
so its cosmological effects are observed solely through its 
influence on the expansion rate, H. Since the expansion 
is accelerating, the dark energy equation of state — 
P<p/p<t> < 0. Assuming that the equation of state varies 
slowly with time, the dark energy component must only 
recently have become cosmologically significant. In order 
to best constrain this component, it is therefore desirable 
to probe the redshift range < z < 2. 

Several authors have suggested using counts of clusters 
of galaxies to probe the evolution of the d ark energy in 
this redshift range (Haiman et al. 2000c iKneissl et al. 



2002; Huterer fc Turner 2001; Podariu fc Ratra 2001; New- 



man et al. 2002). Because these number counts depend 



on the growth function as well as on the proper distance 
probed by Type la Supernovae (SNe) and Cosmic Mi- 
crowave Background (CMB) experiments, a cluster survey 
will yield a constraint complementary to constraints from 
those experiments. Another advantage is that clusters are 
big, bright, and sparse enough to make surveys of large 
volumes relatively tractable. 

Many factors affect our ability to extract constraints on 
cosmological parameters from a cluster survey. In this pa- 
per we examine the trade-offs in parameter estimation be- 
tween conducting a narrow deep survey which may require 
intimate knowledge of the nature of the dark energy, and 
a broad shallow survey which may not. This is because 
shallow surveys use the evolution of the cluster mass func- 
tion to constrain primarily f2 m , and are much less sensitive 
to fl,p or the equation of state, w^. In this case, the in- 
terpretation of the cluster counts will be relatively free of 
the modeling parameters. On the flip side, once several 
of the other parameters are determined reasonably accu- 
rately by combining a shallow survey with CMB and SNe 



constraints, a deep cluster survey will be crucial if one 
hopes to measure any evolution in the equation of state. 

Considering the explosion of cluster survey activity that 
will occur in the next few years, it will be important to ei- 
ther resolve or work around the theory uncertainties in 
order to make the best use of the imminent data. One 
problem that must be examined is how to relate the ob- 
servable quantity (e.g. cluster X-ray gas temperature) to 
the quantity predicted by the cosmological theory (usually 
mass). In this paper we determine quantitatively how un- 
certainties in the Mass- Temperature (M-T) relation and 
in the nature of the dark energy will affect the optimal 
constraints, and in light of that information, identify the 
most effective types of cluster surveys to carry out. Fi- 
nally, accepting the current limitation of the theory, we 
show how the constraints depend on the total number of 
clusters and implicitly on the limiting mass of the survey. 

2. BACKGROUND 

There is a long history of using various cosmological 
tests to constrain parameters in different directions of pa- 
rameter space. Often an observation will not be particu- 
larly useful in constraining an individual parameter, but 
will rather nicely constrain a sum or product of parame- 
ters. The downside is that usually there is another com- 
bination of parameters which is very poorly constrained; 
a degeneracy orthogonal in parameter space to the best 
constrained combination. To break the degeneracy, we 
can combine the constraint with other cosmological exper- 
iments. The key element in obtaining a narrow region of 
likelihood contour intersection is finding observables that 
each have a different dependence on the cosmological pa- 
rameters. 

Much progress has been made in constraining the rather 
large set of parameters typically present in today's cosmo- 
logical models. One of the most import constraints comes 
from observations of the CMB. Existing CMB experiments 



(Smoot et al. 1992; Miller et al. 1999; de Bernardis ei 
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al. 2 001; |Hanany et al. 2000] |Padin et al. 200l| ; [Pryke 
et al. 2001) have already yielded information about the 
initial conditions at the time of last scattering, and have 
made some progress in determining the physical matter 
and baryon densities: f2 m /i 2 , flhh 2 , and the angular dis- 
tance to the last-scattering surface rg. Now with the up- 
coming MAP and Planck missions, the CMB anisotropics 
will be studied with increasing precision, and the con- 
straints will be improved by incorporating new data on the 
polarization of the CMB. The three primary parameters 
mentioned above will be measured to accuracies approach- 
ing 1%, with auxiliary parameters including the equation 



of sta e w,!,, measured to accuracies a few times worse than 



this ( Tue^T 



hence the matter density f2 m (in units of the critical den- 
sity Pcrit = 3H§/8nG, where H = lOO/i kms" 1 Mpc^ 1 
is the Hubble constant). In models with a lower matter 
density, the growth of structure ceases earlier. When nor- 
malized to the same level of fluctuations today, this implies 
larger fluctuations at earlier times. If the number density 
of objects is related to the amplitude of fluctuations on a 
length scale R oc M 1 / 3 , as in the Press-Schec hter theory 
or its generalizations ( Press fc Schechter 1974 ), then this 
predicts more clusters of a given mass at high-z in models 
with a low f2 m . Th is method has been tested numerically 
by many authors (|Efstathiou et al. 1985 , |Efstatliuou fc 
Rees 1988|, |White et al. 19931, |Lacey fc Cole 1994 IGelb fc 



Bertschingcr 1994L |Bond fc Myers 19961, |Eke et al 1996 



The constraint on rg can be interpreted as essentially 
fixing one combination of the matter and dark energy den- 
sities and the evolution of the dark energy. However, even 
when we impose the condition £lk = 0, there are still de- 
generacies thats exist between J7 m and the Hubble param- 
eter, and to a smaller extent, the equation of state of the 
dark energy, w^. 

There is a well known complementarity between the 
CMB constraints and constraints from type la supernovae 
observations which one can take advantage of to break 
these degeneracies ( White 1998 L Tegmark et al. 1998|, Per 



Shcth fc Tormcn 1999| , and |Jcnkins et al. (200 1)|) , and 



mutte ret al. 1999] , Hu et 



al. 



1999| ). The supernova 
constraints also come from a distance measure, but in 
parameter space they are orthogonal to the constraints 
from the CMB because the distance involved is much 
smaller than the distance to the last scattering surface, 
from which all of the CMB photons originate. Thus for 
example, in measuring the energy densities of the matter 
and the dark energy, the constraint from the Type la Su- 
pernovae measures f2 m — Qd, . whereas t he CMB measures 
a weighted sum of fl m + Q$ ( |Cohn 1998 ). The effect is not 
quite as pronounced in the £l m — w,p plane, but it seems 
clear that supernovae will provide a strong complementary 
measurement. The intersection of the likelihood contours 
from a supernova experiment such as SNAP with CMB 



consti aints obtained using PLANCK could determine the 
matter densit y to 3.3% it the polari zation information is 
included (e.g. Tegmark et al. 1998). An analysis of the 
optimal survey strategy for Supernovae searc hes has been 
carried out by Spergel & Starkman (2002), who deter- 
mined that intermediate redshifts of 0.3 — 0.4 are the best 
ones to probe when combining the Supernova results with 
the CMB. 

Although supernovae do seem to provide a clear way to 
narrow the parameter space, it would be a very valuable 
cross check if a third observable was found whose likeli- 
hood contours intersected the other two in the same re- 
gion of parameter space. Such an independent constraint 
could additionally remove one of the largest degeneracies 
in determining the equation of state and its time evolu- 
tion. A promising observable to study is the number den- 
sity of rich galaxy clusters as a function of rcdshift, an idea 



which was su ggested as early as 1992 by Oukbir & Blan- 



chard (1992) among others. Studying the number density 
would provide an orthogonal constraint because in addi- 
tion to being sensitive to the comoving distance measure, 
it is also sensitive to the rate of growth of fluctuations and 



it was found that the formalism accurately predicts the 
abundance of rich clusters as a function of redshift for a 
broad range of different cosmologies. 

Unfortunately, there are still many factors that present 
difficulties when attempting to use a cluster survey to 
constrain the cosmolo gical parameters. A recent study 
( Viana fc Liddle 1999 ) concluded that when all of the ma- 
jor sources of error both in theory and observation are ac- 
counted for, an unambiguous determination of the matter 
density f2 m is not yet possible. One complication is that 
the model of the evolution of the dark energy equation 
of state is not well determined by any other observations. 
This makes it difficult to interpret the impact of observed 
evolution of the number density with redshift on the the 
parameter Q m , because it becomes entangled with w^(z). 

Upcoming studies will introduce a golden age in cluster 
surveys. Several observational programs will replace the 
current catalog of tens of clusters at intermediate redshifts 
(z < 1) with thousands of detections improving not only 
the quantity, but the quality of the available data. With 
the launch of the XMM/Newton and Chandra satellites, 
the study of X-ray emission from the IGM has dramati- 
cally advanced. The XMM (X- ray Multi Mirror Satellite) 
Large-Scale Structure Survey (Pierre, Valtchanov & Re' 



frcgier 2002j ) will map the location of clusters and groups 



of galaxies out to z w 1 over an 8° x 8° area. This sur- 
vey alone is expected to detect up to 800 clusters in the 
temperature band ranging from 0.4 to 4 keV. Another pro- 
posed survey usi ng the same instrum ent, the XMM Clus- 
ter Survey (XCS; [Romcr et al. 2001 ), is expected to cover 
pa 800 deg 2 and include more than 8000 clusters, as many 
as 1800 of which will have well-defined temperatures (es- 
timates for n m = 0.3, fi = 0.7) with T > 2 keV. There 
should even be a significant number detected with z > 1. 

Meanwhile, surveys for clusters are now being carried 
out at optical wavelengths (using both galaxy counts and 
weak lensing) and using th e Sunyaev-Zel'dovich effect 
(Sunyaev & Zel'dovich 1980). In the more distant future, 
the Planck satellite has the potential to find 10 4 clusters 
(with some dependence on cosmology) out to z ~ 1. This 
enormous quantity of upcoming cluster survey work mo- 
tivates our study to find survey parameters that optimize 
the constraints on cosmological parameters. 
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3. METHOD 

In order to address these q uestions we make use o f the 
Fisher matrix formalism (sec Tcgmark et al. 1996| for a 
pedagogical review) . This allows us to elucidate the depen- 
dence of our hypothetical constraints on both our survey 
design and our theoretical uncertainty, highlighting the 
role of parameter degeneracies. Because we have chosen 
to use the Fisher matrix, expressed in terms of the deriva- 
tives of the likelihood function at the maximum, rather 
than a full likelihood analysis, all of our constraints are 
"local" rather than global. This should not affect any of 
our conclusions. 

To make our analysis more relevant to the future where 
CMB and SNe surveys have progressed significantly, we 
consider the constraints from a combination of hypothet- 
ical observations. First we have a cluster survey, but we 
include as 'prior' information the Fisher matrices obtained 
from a mock MAP mission and a future SNe survey. We 
describe each of these pieces below. 

3.1. Assembling the cluster Fisher matrix. 

We assume a flat cosmology and have chosen as our basis 
of parameters the space (f2 m = 0.3, h = .67, fl^h 2 = 0.018, 
n s = l,w <j> = -1, 5 H = 5.02 xlO" 5 , T/S = 0, r = .2, M-T 
normalization = 6.59 x 10~ 10 , <jm-t — -1). The quan titie s 
M-T normalization and <jm-t are defined below in §3.2. 



If we bin narrowly in the observable temperature and in 
rcdshift, the abundance of clusters is naturally described 
by a P oisson distribution, and so the Fisher matrix F be- 
comes (Holder et al. 2001) 



8 2 \nL 
' 86,89, 



^ 8N a 8N a 1 
^ 89, 89,j N a 

a—1 J 



(i) 



where the sum runs over S bins, chosen to be sufficiently 
small that our fiducial model predicts N a -C 1 clusters per 
bin. In practice, binning in redshift finer than Az = 0.01 
does not change the constraints significantly. Here C is the 
likelihood function, the 9i are the cosmological parameters 
and N a is the number of clusters expected in bin a. The 
derivatives are evaluated at the position of the fiducial 
model in likelihood space logarithmically via the rule: 



9 



x AO) - g(6i A0) 
29, \n(A9) ' 



(2) 



All derivatives were found to converge with A9 = 1.0001. 
When varying Q m we also varied fl\ to maintain Qfc = 
0. The dependence on the cosmological parameters enters 
through th e mass function dn/ dM. We use the fitting 



function of Jenkins et al. (2001) 



dn 
dM 



diner 



M 2 din ill 



0.315 exp 



(-lncr + 0.61) 



3.8 



(3) 

where p m is the mean background matter density, and a 2 
is the variance of the mass, filtered on some scale R(M). It 
contains all the details of the cosmological model through 
its dependence on the power spectrum P(k): 



Airk 2 dk W 2 (k)P(k) 



(4) 



We have chosen a spherical top hat (in real space) as our 
smoothing window, with Fourier transform W. We use a 
power-law primordi al power spectrum and the fits to the 



transfer function of 



Eisenstein fc Hu (1997) 



Using the M-T relation, along with its associated Gaus- 
sian scatter, we map the number densities by mass, n(M), 
to number densities by temperature, n(T). We construct 
the Fisher matrix only from bins with a temperature 
greater than a certain threshold; this corresponds to the 
detection threshold of actual surveys. Though the major- 
ity of our analysis is done with a threshold of 5 keV, we 
also explore the tradeoffs behind choosing different tem- 
perature thresholds. One benefit of this method is that 
since the volume at a given redshift scales linearly with 
angular survey size, so too does the number of clusters ex- 
pected for a given temperature range. Therefore, we have 
the relation N a = n a V, where V represents some num- 
ber of (/i _1 Mpc) 3 and n a is the number of clusters per 
(/i _1 Mpc) 3 . This identity allows us to rewrite (II]): 



F -yV dria dna 1 

a—1 J 



(5) 



As V is proportional to the angular survey area, multi- 
plying the area by some factor corresponds to multiplying 
each Fisher matrix element by the same factor. We will 
take advantage of this below. 

3.2. Constructing the M-T relation 

For the cluster M—T relat ion we take the 'virial' relation 
dKaiser 19861 |Peacock 1999|) 



TocMX 3 H 2 / 3 (z)A 1 / 3 (z)h 



2/3 



(6) 



which can be derived by assuming clusters are self-similar, 
pressure supported, isothermal spheres. Here M v - n repre- 
sents the virial mass of the cluster, A(z) is the mean inte- 
rior overdensity, relative to the critical density, for which 
we say the object has virialized. A is the familiar factor 
187r 2 for an Ct m — 1 univers e, and otherwise is given by 
the fitting function found by Bryan fc Norman (1986)| . 



We fix the normalization of the M-T relation by appeal- 
ing to numerical simulations. However, the normalization 
constant could potentially vary across different cosmolo- 
gies. A comprehensive summary of all the different cal- 
cula tions of the normalization is fou nd in Henry (2000) 
and Pierpaoli, Scott, & White (2001). Since we will com- 



pensate for differences in h, the parameters we should be 
concerned with are f2 m and ilb- We will use the value 



quoted by Eke et al. (1998) (hereafter ENF), 



T = 6.59 



M 



2/3 



keV 



(7) 



with a Gaussian scatter of om-t = 13.33%, because the 
cosmology they simulate is very close to our fiducial model 
(fi m = 0.3, n A = 0.7, h = 0.7, n h = 0.04). Other than the 
difference in h, we match these exactly. This normaliza- 
tion was calculated by ENF as an average X-ray emission 
weighted temperature of 10 clusters at z = 0. 
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Combining the normalization with the proportionality 
wc arrive at the final form of our M-T relation: 



2/3 



(fi m ,o(l + zf + n 0)O (l + zf^) V3 , (8) 

where A(0) denotes the value of A(z) for the fiducial 
model, and is not recalculated when we vary the cosmol- 
ogy to calculate derivatives. For now, we will assume that 
the evolution of the M — T relation is described (perfectly) 
by Eq. (j|) over the range of redshift of relevance. We shall 
return to this point later. 



3.3. Assembling the CMB and SNe Fisher matrices 

To make our analysis more realistic, we have computed 
the Fisher matrices associated with CMB and SNe exper- 
iments, and added them as priors to the Fisher matrix 
obtained with a cluster survey. To compute the Fisher 
matrix associated with an experiment such as the upcom- 
ing MAP mission, we have used the methods described in 
Eiscnstein et al. 199£ and have faithfully reproduced their 



constraints. Assuming Gaussian perturbations and noise 
the CMB Fisher matrix is given by 



X.Y 



1 

XY 



oc 



Yf 



(9) 



Cxi is the power in the I th multipole and the indices X 
and Y run over T,E,B,C; the temperature, E-channel po- 
larization, B-channel polarization and the temperature- E 
mode cross correlation. The matrix Cov^ depends on the 
Cxi, the beam window function, and detector n oise levels. 
The explici t form of the matrix can be found in Eiscnstcir] 
et al. 1999 where we have used the specifications for MAP 



at 90 GHz. We have optimistically assumed a sky coverage 
of 80%. 

Calculating the Fisher matr ix for a SNe experime nt is a 



problem addressed recently in Tegmark et al. 1998 . It de 



pends on several observables; the number of SNe observed, 
N, the mean redshift of the observed SNe, z, the scatter 
about this mean redshift, Az, and the standard deviation 
in the observed magnitude from the predicted value (Am) 2 
. Because we would like to use the constraints that will 
be obtained from future experiments such as SNAP, we 
also postulate a Gaussian distribution function f(z) of N 
supernovae distributed about z. For a flat cosmology, the 
Fisher information matrix is given by 



1 



where 



(Am) 



Wi(z) 



f(z)w i (z n )w j (z n )dz, 



( 5 



-i dy 

d0; 



(10) 



(11) 



The 9i are the cosmological parameters to be constrained, 
and rj is given by 



da 



a a 2 H(a) 



(12) 



where H{a) = H [n nl a- 3 + n <p a^ 1+w ^} 1 ^ . 

We computed the SNe Fisher matrix for two different 
cases. First, we calculated the constraints given by a cau- 
tious estimate of 200 supernovae. They followed a Gaus- 
sian distribution at a mean redshift of 0.65 with a dis- 
persion of 0.3 in redshift. The standard deviation in the 
observed magnitude was Am = 0.3. We also considered a 
more optimistic survey, providing 400 supernovae with a 
mean redshift of 0.7 and a dispersion of 0.4 in redshift and 
Am = 0.2. Both of these scenarios are relatively conserva- 
tive when compared to the expected yield of around 2000 
Supernovae to be observed by the proposed experiment 
SNAP. Both samples are used in the analysis below. 

4. RESULTS 

It is instructive to compare Fisher matrices obtained us- 
ing different experimental setups. In this paper, we change 
the specifications of the cluster survey, and find the effect 
on the constraints obtained. An important question to 
consider is whether the limitations on improving a con- 
straint are dominated by the experiment or the underly- 
ing theory, or some combination thereof. In the analysis, 
we have separated the two issues by parameterizing the 
theoretical uncertainty in the M-T relation, and including 
its normalization and scatter in the cluster Fisher matrix. 
To study experimental effects alone, we simply invert the 
Fisher matrix without any additional priors on the M-T 
relationship. To examine the effect of theory uncertain- 
ties, we add additional priors for the M-T normalization 
and scatter for a given experimental setup. Thus, we gain 
an intuition for which sources of error are the most signif- 
icant. 

An observational question that merits attention when 
designing the optimal cluster survey is the trade off be- 
tween angular coverage on the sky, and depth of field of 
the survey, z G nd- It would be nice to know whether the 
bulk of the constraint on a parameter is determined by 
the redshift evolution of the number density of clusters, or 
rather simply on the total number of clusters today. The 
former would suggest that a narrow, deep survey would be 
best, while the latter could easily be accomplished with a 
shallow survey with a wide angular coverage. In general, 
the constraint depends on both quantities, and the domi- 
nant effect depends upon the parameter being considered. 
We approach the question from the back by simulating 
both types of observations and plotting the dependence of 
the constraint on depth or angular coverage. 

In our analysis of the various cluster surveys, we wish to 
consider not only the survey's constraining capabilities on 
its own, but also its ability to improve the constraints al- 
ready obtained from observations of the CMB and type la 
Supernovae. To obtain combined constraints, one has sim- 
ply to add the Fisher matrices together and invert the re- 
sulting matrix. In p ractice, using Cholesky decomposition 
( Press et al. 1986| ) to calculate the inverse is somewhat 
unstable because the fiducial parameters vary by many 
orders of magnitude. This instability is avoided by per- 
forming the analysis on the logarithms of the cosmological 
parameters. Since the Fisher matrix transforms as a co- 
variant tensor under a variable change in parameter space, 
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this modification is easily accomplished using the Jacobi 
matrix. 

The ith. diagonal element of the inverse of the sum of 
the three Fisher matrices gives a lower bound on the 1 — a 
error bars obtained from all three experiments. When the 
other parameters are allowed to vary: 



A0? > (F- 1 ), 



(13) 



This relation accommodates the propagation of errors in 
different parameters, since inverting the Fisher matrix is 
equivalent to marginalizing over all of the other parame- 
ters. We define the optimal percent constraint as the min- 
imum A9i/0i allowed by equation 13, which is interpreted 
as the size of the 1 — a error bar. 

However, the Fisher matrix also quantifies the covari- 
ance of the parameters. For each type of observation, a 
realistic picture of its capabilities is obtained by plotting 
elliptical curves of constant likelihood using a 2 parameter 
subspace of the inverted Fisher matrix. This view explic- 
itly depicts the degeneracies discussed in §|| between the 
two displayed parameters, while marginalizing over all the 
rest. These degeneracies can be broken if the curves of 
constant likelihood from other experiments are oriented in 
an orthogonal direction, since adding the Fisher matrices 
effectively multiplies the likelihoods. Figure 1 displays the 
l — o" regions of likelihood for all three individual experi- 
ments, and also the combined constraint, illustrating the 
efficacy of this technique. Table 1 contains the \ — a errors 
on each of the parameters after combining the CMB ma- 
trix (using specifications of MAP with polarization), 200 
supernovae matrix (about a mean redshift of 0.65), and a 
cluster survey of 1000 deg 2 with temperature threshold 5 
keV to a limiting redshift of z on d = 0.7. 

Table 1 reveals an interesting phenomenon; the X-ray 
cluster survey constrains the normalization of the M-T 
relationship, if its redshift evolution is assumed known. 
Without any prior information about the normalization of 
the the M-T curve from simulations, a constraint of reason- 
able order can be achieved observationally. The constraint 
on the scatter gm-t is not as good. This indicates that 
unless theory can predict the M-T normalization to better 
than ~ 20%, prior knowledge of the normalization will not 
significantly alter the results (see e.g. Pierpaoli, Scott, & 
White 2001). Prior knowledge of the scatter appears to 
be even less significant. We quantitatively investigate this 
below (see Figure 7). 

The orientation and size of the cluster ellipse in Figure 
[l] depends on the depth and angular extent of the sur- 
vey. The change in orien tation in the Q m — plane has 
been noticed previously (Hu &Kravtsov 2002). That the 
size will decrease with increasing depth is apparent; more 
information will result in a better constraint. Figure |^ 
demonstrates how the orientation of the ellipse changes in 
a progression from a shallow survey to a deep one at fixed 
angular size. These plots do not contain information from 
either the CMB or the SNe. However, the number density 
of clusters at a given redshift is a quantity which is largely 
insensitive to the parameter Vl^h 2 . This causes the Fisher 
Matrix to be almost singular, reflecting the large degener- 
acy in that direction. Thus we do rely on results from Big 



Parameter 


Fiducial value 


1 — ct error 
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\J.O 


lU.O/0 


h 


67 




O, h 2 




O.OU /0 


ns 


1 


Z.JO /0 




-1 


8.90% 


Sh 


5.02 x 10~ 5 


5.05% 


T/S 





± 0.796 


T 


0.2 


17.5% 


M-T norm 


6.59 x 10- 10 


15.4% 




0.1333 


138% 



Table 1 

The errors in each parameter from a combination of a CMB, 
SNe, and cluster survey of 1000 deg 2 to a depth of z=0.7. 
Percent errors are given in terms of the fiducial values. 



Bang Nucleosynthesis to provide us with a 10% prior on 
this quantity, which we include throughout our analysis for 
greater numerical stability. Note that the shape converges 
for the deepest surveys, as fewer and fewer clusters above 
the detection threshold are detected. 

The change in the orientation of the direction of degener- 
acy is a direct consequence of the fact that cluster surveys 
probe a combination of the proper distance and the growth 
function. In the range < z < 0.5, the proper distance 
to adjacent redshift bins is a quickly varying function, and 
thus the sensitivity to small changes in the cosmological 
parameters is dominated by changes in this proper dis- 
tance. Hence, the direction of degeneracy of these shallow 
cluster surveys tend to align with contours of constant 
proper distance, as can be seen in the first panel of Fig- 
ure 2. In the high redshift range z > 1 however, each 
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Fig. 1. — 1 — ct regions of maximum likelihood in the Q m — w 
plane. Constraints come from a cluster survey of 1000 deg 2 to a final 
redshift of 1.2 and temperature threshold 5 keV (solid), observations 
of the CMB (dashed) using the specifications of MAP (outer) and 
PLANCK (inner) including the polarization information, and finally 
a hypothetical supernova observation (dot-dashed) with 200 (outer) 
and 400 (inner) supernovae. The solid contour is the joint constraint 
from the cluster survey, 200 supernovae, and MAP. 
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rcdshift bin adds an increasingly large number of clusters. 
The number of collapsed objects at high redshift is very 
sensitive to the growth function, and thus in this redshift 
regime, the direction of degeneracy aligns instead with the 
contours of constant growth factor, normalized to today. 
This leads to an interesting detail: populations of different 
masses will come to be dominated by the growth factor at 
different redshifts. This is because of the hierarchical na- 
ture of structure formation. More massive clusters form 
later, and thus massive clusters get rarer faster as deeper 
redshifts are probed. The more massive populations get 
dominated by the growth function at smaller redshifts. We 
have determined that the higher the temperature thresh- 
old, the lower the final redshift z on( j at which one runs com- 
pletely out of clusters; z C nd ~ 1 for 8 keV and z cn d ~ 1.6 
for 5 keV (see Figure 2), while z cn d is far beyond 2 if one 
stretches the threshold to 3 KeV. The orientations of the 
constraints at these limits are also slightly different, being 
most orthogonal to CMB and supernova constraints for 
the smaller objects. One complication is that the lower 
the temperature threshold is extended, the more one has 
to worry about extra physics that may affect the low mass 
systems, but are difficult to theoretically model. These 
considerations must be traded off with depth and angu- 
lar coverage when a survey is being planned. The optimal 
strategies clearly depend on the capabilities of the survey- 
ing technique, which we discuss below. 

By far the easiest experimental specification to examine 
is the angular coverage of the survey. Changing the an- 
gular coverage effectively changes the volume element at 
each redshift z. As the Fisher matrix contribution from the 
cluster survey is linearly dependent on its volume, and the 
volume is, in turn, linearly dependent on the survey's an- 
gular size, the parameter constraints solely from the clus- 
ter survey go as: 

M K 7§ (14) 

where S is the solid angle subtended by the survey, where 
for a full sky survey S is 47r/41000. The volume element 
also depends on the comoving distance, which is a func- 
tion of the cosmological parameters. Figure 3 shows the 
improvement in the constraint on the matter density fi m 
for the three experiments as the angular size of the survey 
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Fig. 2. — 1 — a regions of maximum likelihood in the Q m — 
Wfj, plane for a series of cluster surveys differing only in survey 
depth. In the left hand panel, constraints are plotted for z en d = 
0.1, 0.2, 0.3, 0.5, 0.7, 0.9 decreasing in size. The right hand panel has 
z end = 0.9, 1.1, 1.3, 1.5, 1.7, 1.9. 



increases. The relation is not quite as simple as in equa- 
tion 14 because of the prior information provided by the 
CMB and Supernovae. 

What we notice is that if the survey is deep enough there 
seems to be a "sweet spot" , a clear cutoff place where the 
effort involved in performing the survey returns maximum 
scientific information, and beyond which the payoff is not 
high. For example, for survey depth z on( j = 1.2 the payoff 
for extending the observations from 500 to 2000 deg 2 is 
quite considerable, but extending the survey to 10 4 deg 2 
does not improve the optimal constraint significantly from 
a survey of <~ 2000 deg 2 . It is important to keep this 
turning point in the back of one's mind because increasing 
the angular coverage greatly increases the total number of 
clusters that need to be observed. 

Another strategy to consider is to make the survey 
deeper rather than wider. In this case there may be a 
smaller total sample of clusters, but the evolution of the 
number density will be more apparent. In Figure 4, the 
dependence of the constraint on the matter density is dis- 
played as a function of z cn( j, the depth of the survey. 
Of course many surveys are flux limited, and thus the 
depth of the survey is different for different mass scales; 
brighter more massive objects can be seen at greater dis- 
tance. However, since a significant portion of the con- 
straint comes from the evolution of the density of clusters, 
and not just from the total number, it is crucial that the 
sample be complete. Therefore, the depth of the survey 
is effectively the one at which completeness is achieved 
for the faintest objects being counted, and we will assume 
this in our analysis. In practice this depth will depend 
on the type of survey, and of course on the temperature 
threshold T m i„. Constraints in Figure 4 are shown from 
the cluster survey alone, the cluster survey combined with 
the expected constraint from the CMB MAP mission (with 
polarization information) , and finally from the cluster sur- 
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Fig. 3. — The percent constraint on Qm as a function of the 
angular coverage of the cluster survey. We have added the fisher 
matrices from the MAP CMB experiment, and from 200 supernovae, 
as well as a 10% prior on Q. h h? from BBN. The solid line is a survey 
to depth 2 en( j = 1.2, the dashed line is a survey to z en( j = 0.7, and 
the dot-dashed line has a depth of z en( j = 0.3. 
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vey combined with both the CMB and SNe experiments. 

Again what we see is that there is a definite "sweet 
spot." The role of the prior information matrices of the 
SNe and CMB experiments is also apparent in this graph. 
The flat nature of both of the dot-dashed curves between 
z C nd — and z cn( j ~ 0.5, which incorporate the constraints 
from all three types of experiments, are a result of the fact 
that although the survey is getting deeper and the num- 
ber of clusters detected is growing, the constraint from the 
cluster survey is weak compared to combined constraints 
from the CMB and SNe. The solid curve displaying the 
cluster constraint alone shows that this is true. In terms 
of the elliptical plots above, the likelihood contours from 
CMB experiments would be entirely within the region of 
maximum likelihood determined by the cluster survey. Af- 
ter z en( j = 0.5 we see a marked improvement, which means 
the likelihood contour will have become small enough to 
intersect the CMB and SNe contour in a way that de- 
creases the most probable region compatible with all the 
experiments. Not only has the size of the cluster constraint 
decreased significantly by z cn( j ~ 0.5, but as one can see 
from Figure 2 the orientation has passed the vertical and 
becomes increasingly orthogonal to the other constraints 
between 0.5 < z cn d < 1.7. The role of the SNe constraint 
is also clear. 400 supernovae provide a narrower lever arm 
in parameter space than 200. Since both SNe contours are 
at a significant angle to the cluster constraints, and since 
they intersect at the position of the fiducial model (by 
fiat), then the effect is simply to decrease the area of the 
intersecting region by an approximately constant amount, 
which can be observed in Figure 4. 

Another interesting property of experiments that should 
be carefully studied is the effect of changing the detection 
threshold temperature on the quality of the constraints 
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Fig. 4. — The percent constraint on Q m as a function of the 
depth of the cluster survey. This survey has an angular coverage of 
1000 deg 2 . The solid line is the constraint from the cluster survey 
combined with the 10% BBN prior. The dashed line is obtained 
when the results of the CMB (MAP), cluster observations, and BBN 
prior are added, and finally the dot-dashed lines are the constraints 
obtained when all three types of observations are used; the cluster 
survey, the CMB, and either 200 (upper) or 400 (lower) SNe, also 
with the BBN prior. 



that can be obtained. In some types of surveys this could 
be intimately related to the depth of the survey, for ex- 
ample if the survey was flux limited, since flux correlates 
with temperature. However this is not always true. In the 
case of the SZE effect, the limiting temperature is largely 
independent of redshift. Also for X-ray surveys one could 
argue that those clusters for which accurate temperatures 
can be known are well above the flux limit, allowing es- 
sentially volume limited samples to be constructed. Since 
temperature can be used as a proxy for mass, studying 
this dependence is akin to asking whether most of the con- 
straint comes from the rare high mass clusters, or whether 
the large number of lower mass clusters are needed in or- 
der to adequately constrain the parameters. We would 
prefer the former scenario, because the higher T m i n , the 
less sensitive are the results to non-gravitational physics 
in the modeling, and the easier it is to make the observa- 
tions. The value of T m ; n has a strong nonlinear effect on 
the number of objects that will be found within the do- 
main of the survey, as Figure 5 shows. Figure 6 illustrates 
the improvement in our knowledge of the matter density 
as we lower the threshold temperature. 

This clearly demonstrates that it is indeed the rarest 
high mass clusters that contribute the bulk of the con- 
straint. For a survey out to z cn< \ = 1.2, if one drops the 
threshold halfway down the range of the plot, from 8 keV 
to 5 keV, then one adds about 800 clusters and improves 
the constraint by a factor of 3/4 the original value. If one 
considers the bottom half from 5 keV to 2 keV, one adds 
~ 24, 000 clusters, yet only improves the constraint by an- 
other factor of 3/4. If all clusters were equally important 
in constraining Q m then the constraint would drop signif- 
icantly more in the bottom half of the plot, because we 
have increased the number of clusters surveyed by more 
than an order of magnitude. 

It is useful to probe how strongly the constraint on Q m 
(or any other parameter) depends on the theory implicit in 
the survey. Many cluster surveys involve observation of the 
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Fig. 5. — The total number of clusters detected at all redshifts 
out to z cnt } = 1.2 (solid), z on d = 0.7 (dashed) and z cn d = 0.3 (dot- 
dashed) as a function of the detection threshold in keV for a fixed 
angular aperture of 1000 deg 2 . 
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X-ray luminosity or temperature, either because they are 
X-ray surveys, or else they are optical or SZE surveys that 
require an X-ray followup to obtain redshift and temper- 
ature information. For a cluster survey in X-ray tempera- 
ture, this uncertainty is framed in terms of the uncertainty 
in the normalization and the scatter about the mass tem- 
perature relationship. It is conceivable that adding more 
clusters to the survey will not affect the constraint on f2 m 
simply because the Fisher matrix is dominated by the er- 
ror from the M-T normalization and scatter. If this is the 
case, considerable effort should be invested in improving 
these theoretical considerations, so as to use the plethora 
of upcoming data to maximum advantage. 

By including the M-T normalization and scatter as pa- 
rameters in the Fisher matrix, we have investigated this 
important issue. We find that precise knowledge of the 
scatter about the M-T relationship does not significantly 
improve the quality of the constraint on Q m . Reducing the 
prior on the scatter in a survey to z on( j = 0.7 from a very 
loose prior of 100% to an extremely tight 5% improved 
the constraint on the matter density less than 1%, with a 
10% BBN prior and the CMB and 200 SNe experiments 
included. The improvement of the constraint on fi m with 
better knowledge of the M-T normalization is of the same 
order when the prior is changes from 10% to 1%. Weaker 
priors on the normalization have little effect because the 
cluster survey constrains it to 15% by itself (Table 1). 
Figure 7 shows this dependence for several redshifts. 

5. CONCLUSIONS 

Many exciting surveys of galaxy clusters have been pro- 
posed, involving several different observational techniques 
to locate the clusters. In this paper we have analyzed their 
potential to constrain two of the most fundamental cosmo- 
logical parameters, f2 m and w^, within a larger framework 
of information gathered from other types of observations. 
We have considered the effects of theoretical uncertain- 
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Fig. 6. — The percent constraint on Q m as a function of the 
temperature detection threshold. We have added the same priors 
as in the previous two figures. The solid line is a survey to depth 
2 cnd = 1-2, the dashed line is a survey to z en( j = 0.7, and the dot- 
dashed line has a depth of z en d = 0.3. 
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Fig. 7. — The percent constraint obtained for Q m as a function 
of the prior knowledge of the error bar on the normalization of the 
M-T relation, also reported in percent. The solid line is a survey 
to depth z cn( j = 1.2, the dashed line is a survey to z en d = 0.7, and 
the dot-dashed line has a depth of z cn( j = 0.3. We have also shown 
the change of constraint as a function of the prior knowledge of the 
scatter from the M-T relation for z cn d = 1.2. The constraint does 
not depend strongly on this quantity. 



ties on these constraints. Even when the equation of state 
Wrf, is postulated to be a constant, large degeneracies ex- 
ist in the f2 m — subspace from all types of observations 
studied here. We have varied some of the characteristic at- 
tributes of cluster surveys such as the limiting temperature 
threshold, the depth of the survey, and the angular cover- 
age, and have noted that the direction of elongation of the 
elliptical region of maximum likelihood depends consider- 
ably on these attributes. Since we know the direction of 
degeneracy for other important experiments, we can use 
this information to optimize the cluster survey strategy, 
taking into account the various limitations for X-ray fol- 
lowup observations. 

We have deduced that in order to achieve maximal sen- 
sitivity to r2 m , with minimal sensitivity to w$ (without 
taking into account CMB and SNe experiments), a fairly 
shallow survey to 0.4 < z cn d < 0.7 (depending on the limit- 
ing temperature used) with large angular coverage is an ef- 
fective approach, because in this case the constraint ellipse 
is vertically oriented in the f2 m — plane. It would be in- 
teresting for future work to see if such a constraint would 
be relatively independent of wL. If so, such a measure- 
ment would yield a significant lever arm to experiments 
such as SNAP, who would then be in a good position to 
tighten constraints on wL. We also found that for a sur- 
vey of fixed angular extent, the clusters with the highest 
constraining power are th ose with high temp erature. We 
confirm the claim made by Holder et al. 2001 that it is the 
clusters past z > 0.5 that provide most of the constraint 
on f2 m (given our assumption that Q m + = 1). This 



is especially evident when SNe and CMB experiments are 
combined with the cluster survey, because the part of the 
cluster constraint that is orthogonal to the SNe and CMB 
ellipses comes from the clusters located at z > 0.5. We 
also see a natural cutoff in limiting redshift beyond which 
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the constraint on f2 m stops improving, as the number of 
clusters hotter than the detection threshold falls off steeply 
with z. More surprisingly, we also find that when the clus- 
ter survey is combined with the CMB and SNe constraints, 
increasing the field of view of the survey does not improve 
the results indefinitely. 

In summary, when designing an experiment, several 
competing effects must be accounted for. If one is in- 
terested in constraining Q m and w^, simultaneously, then 
maximizing orthogonality with the CMB and Supernova 
constraints is desirable because it will maximize the re- 
turn for each unit of observing time that is spent counting 
objects. In general this suggests that deep surveys are ad- 
vantageous because it is at high redshifts that the growth 
function significantly affects the number density of clus- 
ters. If an observable such as the SZ effect can truly probe 
out to very high redshifts with relative completeness, then 
it is also advantageous to stretch the temperature thresh- 
old as low as possible, since this allows a deeper effective 
depth before the clusters run out. The resulting difference 
in angle is a relatively small effect however, and naturally 
there are limits to the practicality of lowering the thresh- 
old. The number of dim distant objects that would need 
to be followed up grows exponentially as the threshold 
decreases. To some extent, this can be compensated by 
making the survey narrower. 

If, on the other hand, the survey is flux limited and can 
probe completely only to a relatively shallow z C nd, such as 
an optical or X-ray survey of clusters, the optimal strat- 
egy changes. In this case we take advantage of the fact 
that very massive clusters are dominated by the growth 
function at smaller redshifts. For maximum orthogonal- 
ity in shallow surveys, one raises the temperature thresh- 
old and studies only the most massive clusters. We have 
shown that this will not significantly weaken the survey's 
performance since typically the more massive clusters con- 
tribute the bulk of the constraining power. Of course to 
compensate for the fact that huge cluster are rather rare, 
the survey can be widened, keeping in mind that for each 
enc i, there is a point of diminishing returns when doing 
this. 



In i nvestigating different survey techniques' abilities to 
cffectwely whittle Hnwn the likely region nhtainrH from 



SNe and CMB experiments, we have been conservative 
about the theoretical uncertainties involved in connecting 
the observable (X-ray temperature) to the quantity pre- 
dicted by the theory (cluster mass) . We have included the 
normalization and scatter as parameters in our Fisher Ma- 



trix ahalysis. Any error in these parameters propagates to 
the constraints on the cosmological parameters. We have 
invest igated the sensitivity of our constraints on 17 m and 



w,/, to errors in the normalization and scatter of the M- 
T relation by including prior constraints on these quanti- 
ties and observing the improvement on the parameter con- 
straints. We found that in general, the constraints on r2 m 
and Wfj, are fairly robust against large errors in the scatter 
of the M-T relation, but that they do depend somewhat 
more strongly on errors in the normalization. It is inter- 
esting that if no priors are added, the cluster survey itself 
will constrain the normalization to around 10%, suggesting 
that cluster surveys could lead to some insight into intrin- 



sic cluster properties. We have not addressed what might 
happen if non-Gaussian effects were to account for a large 
portion of the error on M-T, nor have we accommodated 
a normalization that changes with redshift. These would 
change the shape of the temperature-redshift distribution 
and could impact the determination of the cosmological 
parameters profoundly. To asses the impact of the first 
it would be necessary to employ a monte-carlo method, 
rather than the simple Fisher matrix techniques used in 
this paper. To understand quantitatively the question of 
evolution is beyond the scope of this work. 

It is clear that surveys of galaxy clusters will play a key 
role in the future of further constraining the cosmological 
parameters. Tight constraints on £l m and the evolution 
of the dark energy are essential to developing our mod- 
els of the universe. A precise determination of f2 m will 
help quantify the dark matter problem, and aid in solv- 
ing questions of structure formation. The more f2 m can 
be independently constrained, the better leverage we will 
have later when observations of very deep sources begin to 
give us more information about the evolution of the dark 
energy. 
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